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Abstract 

The9>observation  of  cavitation  in  C.S.M.  's  at  the  throat  of  a  Venturi  tube 
provides  an  ideally  simple  and  direct  method  to  monitor  the  liquid  quality  in 

terms  of  the  concentration  of  unstable  nuclei  as  a  function  of  the  applied 

a 

tension.  In  the  previous  work, of  Oldenziel  and  Lecoffre,  the  use  of  C.S.M. 's 
was  limited  to  the  measurement  of  the  tensile  strength  of  the  weakest  nuclei 

rn »  ; 

present  in  the  liquid.  The  purpose  of  the  currefit, research  is  tq» extend?this 
application  to  determine  the  concentration  of  unstable  cavitation  nuclei  as  a 
function  of  the  applied  tension  over  a  wide  range  so  that  the  nuclei  number 
distribution  can  be  derived.  Various  approximations  are  made  to  construct  a 

J 

simplified  but  still  representative  physical  siodel  of  C.S.M. 'sj  these  factors 
included  the-  liquid  sample  initial  conditions  and  properties,  critical  tension 
and  concentration  of  cavitation  nuclei  (i.e.  liquid  quality)  and  required 
accuracy  and  range  of  the  nuclei  concentration  measurement.  The  overall 
performance  of  C.S.M. 's  is  described  in  terms  of  a  set  of  dependent  parameters 
including  the  volume  flux,  velocity  and  pressure  at  the  throat  and  at  the  exit 
(these  are  fluid  mechanical  parameters);  the  throat  volume,  diameter  and 
length,  and  exit  diameter  (geometrical  parameters);  cavitation  rate  and  bubble 
response  time  (cavitation  parameters);  error  due  to  the  interference  of  more 
than  one  nucleus  simultaneously  present  in  the  cavitation  region  and  error  due 


) 
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to  the  finite  time  duration  of  the  observation  (statistical  parameters); 
nuclei  concentration  measurement,  sample  volume  and  sampling  time  (measuring 
process  parameters  in  general).  The  relative  dependence  of  these  parameters 
is  investigated.  A  simple  quasi  one-dimensional  model  is  used  to  analyze  the 
flow.  Viscous  effects  in  the  venturi  tube  are  estimated  by  matching  an 
axisymmetric  boundary  layer  with  an  ideal  flow  solution  in  order  to 
demonstrate  the  presence  of  a  potential  core  in  the  throat  section  where  cavi¬ 
tation  is  to  occur.  The  dynamics  of  a  spherical  bubble  growing  and  travelling 
in  a  semi-infinite  pipe  are  used  to  deduce  both  the  critical  tension  and  the 
response  time  of  cavitation  nuclei.  The  occurrence  of  nuclei  in  the  cavita¬ 
tion  region  and  the  observation  of  cavitation  events  are  treated  as  stochastic 
processes  of  Poisson's  type  in  space  and  time,  respectively.  On  the  basis  of 
the  assumed  nuclei  number  distribution,  the  expected  errors  due  to  nuclei 
interference  and  to  finite  observation  time  are  thus  related  to  the  volume  of 
the  cavitation  region  and  to  the  duration  of  the  observation  itself. 

Finally,  it  is  shown  that  these  parameters  can  be  organized  into  a  func¬ 
tional  network,  leading  to  a  systematic  representation  of  the  C.S.M.  design 
problem. 

The  most  important  deductions  are  that; 

(1)  the  pressure  in  the  cavitation  region  is  an  extremely  sensitive 
parameter  and  cannot  be  directly  measured  with  the  necessary  accuracy. 

It  is  therefore  proposed  to  deduce  the  pressure  from  the  potential  core 
velocity  at  the  throat  obtained  by  using  a  L.D.V. ; 

(2)  the  volume  of  the  cavitation  region  must  be  reduced  as  much  as 
possible  in  order  to  increase  the  applied  tension  without  incurring  in 
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-  large  errors  due  to  nuclei  interference; 

(3)  the  time  response  of  cavitation  nuclei  imposes  a  lover  bound  for 
the  length  of  the  cavitation  region. 

/ 

Preliminary  tests  of  a  prototype  venturi  made  of  acrylic  plastic  shoved 
that  the  flov  in  the  inlet  can  sustain,  vithout  separation,  the  very  steep  and 
large  contraction  suggested  by  design  considerations,  but  that  it  is  extremely 
sensitive  to  separation  in  the  diffuser  section  folloving  the  test  volume. 

INTRODUCTION 

Early  studies  have  shovn  [5]  that  the  maximum  tensile  strength  that 
liquids  can  sustain  is  considerably  reduced  by  the  presence  of  veak  spots, 
generically  called  "nuclei".  Although  their  nature  has  not  yet  been  fully 
understood,  it  has  been  recognized  that  small  gas  bubbles,  certain  types  of 
solid  impurities,  and  gas  pockets  contained  inside  crevices  and  cracks  of  a 
liquid-solid  interface  are  nuclei,  since  they  act  as  preferred  points  for  the 
onset  of  liquid  rupture. 

The  presence  of  nuclei  is  considerably  important  in  many  flovs  of  techni¬ 
cal  interest,  vhere  they  are  responsible  for  originating  or  variously  modify¬ 
ing  the  occurrence  of  cavitation,  a  generally  undesired  phenomenon  vhich 
severely  limits  the  performances  of  hydraulic  operating  machinery.  In  this 
respect  the  number  of  nuclei  vhich  become  unstable  at  a  given  level  of  tension 
is  therefore  a  fundamental  parameter  in  the  attempt  both  to  predict  the  onset 
of  cavitation  and  to  deduce  scaling  lavs  capable  of  extending  the  results 
obtained  from  model  tests  to  full-scale  operation. 
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The  effect  of  the  presence  of  nuclei  in  liquids  is  not  limited  to  cavita¬ 
tion.  Inhomogeneities  and,  in  particular,  micro— bubbles,  considerably  change 
the  optical  and  acoustical  properties  of  liquids  in  a  way  which  also  depends 
in  general  on  their  concentration  and  size  distribution  [6]  .  Such  phenomena, 
*s  well  as  cavitation,  are  of  great  practical  importance  in  many  technological 
applications  of  underwater  acoustics  and  optics. 

Therefore  significant  efforts  have  been  made  in  developing  methods  to 
detect  nuclei  contained  within  liquids,  whose  effects  usually  dominate  the 
ones  due  to  the  nuclei  localized  at  the  flow  boundaries. 

Some  methods,  such  as  photography  or  holography,  which  produce  a  recorded 
image  of  a  sample  volume  in  the  liquid,  allow  one  to  detect  nuclei  whose  size 
is  larger  than  about  10  pm.  They  also  can  discriminate  between  bubbles  and 
particles,  but  have  the  disadvantage  that  results  are  not  immediately  and 
continuously  available  because  the  analysis  of  the  obtained  data  requires  a 
considerable  time. 

Other  methods,  based  upon  optical  or  acoustical  radiation  scattered  by 
the  nuclei,  can  operate  continuously  providing  results  which  are  readily 
available,  but  do  not  allow  one  to  distinguish  readily  between  bubbles  and 
particles.  For  satisfactory  operation,  a  careful  calibration  using  a  liquid 
sample  of  known  characteristics  is  needed.  Non  uniform  illumination  and  finite 
size  of  the  scattering  volume  may  also  generate  systematic  errors  which  are 
fo  take  into  account.  Besides,  when  particles  are  prevalent, 
considerable  uncertainty  is  introduced  by  the  fact  that  both  their  shape  and 
index  of  refraction  are  unknown. 
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Bubbles  in  liquids  can  also  be  detected  from  the  large  attenuation 
experienced  by  acoustical  radiation  which  resonates  with  the  bubble  natural 
frequency  of  volume  oscillation.  Detection  systems  based  on  this  concept  do 
not  require  long  times  for  data  analysis,  but  are  unable  to  detect  particles 
and  need  comparatively  large  and  homogeneous  sample  volumes. 

For  the  purpose  of  cavitation  research,  a  common  disadvantage  of  all  the 
above  methods  is  that  only  size  distribution  and  concentration  of  the  nuclei 
are  measured.  The  value  of  the  critical  tension  that  makes  each  nucleus 
unstable,  as  required  for  cavitation  studies,  must  therefore  be  deduced 
indirectly  from  size  measurements,  which  can  be  done  only  in  the  case  of 
micro-bubbles.  Solid  particles,  when  detected  and  recognized,  can  be  included 
or  not  in  the  measurement  of  nuclei  concentration,  but  no  information  is  pro¬ 
vided  about  their  critical  tension,  yet  it  is  known  that  only  part  of  the 
solid  particulates  are  active  as  cavitation  nuclei. 

To  overcome  these  difficulties,  cavitation  events  occurring  in  a  known 
pressure  field  can  be  detected  and  counted.  Then  applied  to  bodies  of  stand¬ 
ardized  form  and  size,  this  method  has  indeed  the  advantage  of  providing  a 
direct  measurement.  However,  since  many  nuclei  are  simultaneously  exposed  to 
the  tension,  the  response  of  the  weakest  ones  inevitably  dominates  the  flow, 
so  that  only  cavitation  susceptibility  evaluations  can  be  obtained. 

Recently  Oldenziel,  [1],  [2],  [3],  [4],  proposed  a  new  instrument,  the 
Cavitation  Susceptibility  Meter  (C.S.M.),  in  which  the  internal  flow  of  a 
liquid  through  a  glass  venturi  is  used  to  induce  cavitation  at  the  throat  and 
bubbles  are  detected  optically.  The  C.S.M.  developed  by  Lecoffre  at  Neyrtec, 
[13],  is  based  on  the  same  principle,  but  it  utilizes  a  stainless  steel 
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venturi ,  where  cavitation  bubbles  are  counted  by  recording  the  noise  generated 
by  their  collapse  in  the  diffuser  downstream  the  throat  section. 

Various  methods  to  monitor  the  cavitation  nuclei  population  of  water  sam¬ 
ples,  including  the  above  two  C.S.M. 's,  have  been  recently  compared  together 
at  Delft  Hydraulics  Laboratory  [4].  The  results  confirm  that  C.S.M. 's  can 
reliably  measure  the  critical  tension  of  the  weakest  nuclei  and  their  concen¬ 
tration  up  to  a  maximum  of  the  order  of  10  cm-3,  after  which  the  venturi 
throat  becomes  saturated  with  bubbles  and  the  flow  chokes.  Reportedly  the 
Neyrtec  system  can  also  measure  the  nuclei  critical  pressure  in  the  range 
0  bar  to  -2  bar  without  incurring  in  saturation,  provided  that  the  concentra¬ 
tion  of  unstable  nuclei  does  not  exceed  about  1  cm  3  .  Nevertheless  such  low 
values  of  the  pressure  correspond  to  bubbles  having  an  equilibrium  radius 
between  .1  and  1  jim,  whose  concentration  in  most  water  samples  is  expected  to 
be  considerably  larger  than  1  cm-3.  Therefore  it  appears  that,  during  opera¬ 
tion  in  these  conditions,  several  unstable  nuclei  are  simultaneously  present 
in  the  cavitation  region  of  the  C.S.M.  When  this  happens,  the  pressure  in  the 
throat  increases  significantly  with  respect  to  its  nominal  steady  state  non- 
cavitating  value.  The  C.S.M.  is  then  partially  saturated  and  the  actual  range 
over  which  reliable  measurements  of  both  the  nuclei  critical  tension  and 
concentration  are  possible  is  significantly  reduced. 

Yet,  when  compared  to  other  techniques,  C.S.M. 's  seem  to  have  several 
attractive  features: 

(i)  they  provide  a  direct  measure  of  both  the  nuclei  critical  tension  and 
concentration,  thus  eliminating  the  problem  of  the  uncertain  behaviour  of 
particles  as  cavitation  nuclei; 
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(ii)  data  analysis  only  requires  a  comparatively  short  time,  so  that  the 
measure  can  be  carried  out  in  an  almost  continuous  way; 

(iii)  there  is  no  limitation  on  the  size  of  the  smallest  detectable 
nuclei . 

On  the  other  hand,  successful  operation  of  C.S.M. 's  is  severely  limited  by  the 
nuclei  concentration  compared  to  the  other  techniques  mentioned.  But,  because 
the  C.S.M.  is  such  a  promising  device  in  cavitation  research,  we  have 
concerned  ourselves  in  the  present  paper  with  all  of  the  various  considera¬ 
tions  necessary  for  a  reliable  measurement  to  be  made.  Ve  present  these 
thoughts  here  in  the  hope  that  they  will  be  useful  since  there  is  little 
discussion  of  these  considerations  in  the  literature. 

C.S.M.  OPERATION  ANALYSIS 
General  Approach 

The  configuration  of  the  axi symmetric  C. S.M.  venturi  pipe  considered  in 
the  present  study  is  shown  in  Fig.  1.  The  flow  originates  from  the  far  field 
on  the  left  (subscript  o),  representing  the  sampled  liquid,  to  the  exhaust 
region  on  the  right  (subscript  ex).  Both  the  upstream  and  downstream  condi¬ 
tions  are  supposed  to  be  time-independent.  A  volume  flux  Q(t)  of  liquid  flows 
with  axial  velocity  u(x,r,t)  and  radial  velocity  v(x,r,t)  through  the  duct  of 
local  radius  B(x)  and  cross-sectional  area  A(x)  =  n  B2(x),  divided  in  five 
segments;  a  cylindrical  inlet  (subscript  i),  a  contraction  (subscript  c),  a 
cylindrical  throat  (subscript  t),  a  diffuser  inlet  (subscript  di)  and  a  coni¬ 
cal  diffuser  (subscript  de )  .  The  length  L  of  each  segment  of  the  duct  is 
shown  in  the  same  figure;  contraction  area  ratio  is:  Cr  =  At  /  Ai  and  the 
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diffuser  area  ratio  (expansion  ratio)  is:  C  =  A  /  A* . 

The  main  parameters  used  to  describe  the  overall  performance  of  C.S.M. 's 
are  shown  in  the  block  diagram  of  Fig.  2.  The  existence  of  functional  rela¬ 
tions  between  two  parameters  is  indicated  by  arrows  connecting  the  two 
corresponding  blocks.  The  nature  of  these  relations  will  be  discussed  in 
detail  later. 

Some  of  the  parameters  in  the  diagram  represent  given  data,  some  can  be 
considered  as  independent  requirements  and  some  as  resulting  or  dependent  fac¬ 
tors.  However,  the  subdivision  in  such  categories  is  somewhat  arbitrary, 
since  it  depends  on  the  view  point  chosen  to  approach  the  whole  problem  of 
C.S.M.  design.  The  one  implicitly  adopted  in  most  of  the  present  analysis  can 
be  summarized  as  follows: 

-first  it  is  assumed  that  the  sample  initial  conditions  and  properties 
are  given; 

-second  a  nuclei  number  distribution  (i.e.  the  liquid  quality)  of  the 
sample  is  postulated,  which  is  generally  representative  of  typical  water  sam¬ 
ples  of  technological  interest; 

—  finally  a  set  of  conditions  on  the  remaining  dependent  parameters  is 
deduced  for  the  measurement  of  unstable  nuclei  concentration  to  be  possible 
over  the  assigned  or  assumed  nuclei  radius  range  with  the  required  precision. 

This  point  of  view  is  particularly  useful  for  design  orientation.  On  the 
other  hand,  for  performance  assessment  of  a  given  design,  the  obtainable 
precision  in  the  measurement  of  unstable  nuclei  concentration  over  a  certain 
radius  range  is  the  result  of  the  specific  C.S.M.  geometry  which  can  be 
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analyzed  with  the  methods  contained  herein. 


Liquid  Quality  and  Critical  Tension 

In  cavitation  research  literature  the  liquid  quality  of  a  sample  is  usu¬ 
ally  expressed  by  means  of  the  nuclei  number  distribution:  N(Rq)  =  -  dn/dRQ, 
where  n(RQ)  is  the  number  concentration  of  nuclei  whose  equilibrium  radius  is 
not  smaller  than  Rq.  The  nuclei  number  distributions  of  water  samples  of 
technological  interest  spread  in  practice  over  a  very  large  range.  Here,  for 

simplicity,  we  assume:  N(R  )  =  K..  /  R3  ,  where  K . .  =  10~4  to  10-3  m_1  (typi- 

o  dis  '  o  dis 

—3  —1 

cally  =  10  m  )  ,  which  is  in  the  average  representative  of  the  data 

reported  in  the  literature  in  the  radius  range  10  to  100  |im.  Thus  it  follows 
that 


n(Ro) 


(1) 


In  order  to  relate  the  nuclei  number  distribution  to  the  corresponding 
critical  pressure,  we  make  the  hypothesis  that  all  the  nuclei  consist  of 
equilibrium  gas  bubbles,  for  which  the  critical  tension  with  respect  to 
dynamic  stability,  [5],  is  expressed  by: 


1  + 


Po  ~  Py 

2  S  /  R 
'  o 


(2) 


where  pQ  and  RQ  are  the  equilibrium  pressure  and  radius  of  the  bubble,  py  is 
the  vapor  pressure  and  S  the  surface  tension  (see  Fig.  3).  The  corresponding 
"critical "  velocity  at  the  venturi  throat  deduced  from  ideal,  incompressible, 
steady  flow  calculations  is  shown  is  Fig.  4.  Similarly,  the  exhaust  pressure 
coefficient  =  (pex~pv) /(pQ-pv)  is  plotted  in  Fig.  5  for  the  case  pQ  =  1 
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bar  and  various  expansion  ratios.  When  the  equilibrium  pressure  p  increases, 
the  cpex  curves  shift  to  the  left  in  the  diagram. 

Note  that  the  critical  tensions  of  bubbles  in  the  size  range  of  interest 
(say  10  to  100  pm)  are  quite  small.  Also  note,  on  the  other  hand,  that,  on 
the  same  range,  the  radius  of  marginally  stable  bubbles  (critical  radius)  is  a 
very  sensitive  function  of  the  applied  tension.  These  facts  have  several 
important  (and  unfortunate)  consequences  for  the  design  of  C.S.M. 's.  The 
first  is  that  it  is  virtually  impossible  to  measure  these  small  pressures  at 
the  throat  because  of  inadequate  sensitivity  and,  if  done  intrusively,  separa- 
tion  or  cavitation  may  be  induced.  The  second  is  that  the  throat  pressure 
must  be  evaluated  accurately  to  avoid  large  errors  in  the  determination  of  the 
nuclei  number  distribution.  A  third  consequence  has  to  do  with  the  problem  of 
regulating  the  flow  by  varying  the  exhaust  pressure.  The  curves  of  Fig.  5 
deduced  for  PQ  »  1  bar  show  that  very  low  sensitivity  is  obtained,  unless 
diffusers  with  a  very  small  expansion  ratio  are  used.  Besides,  due  to  the 
shift  of  the  curves  to  the  left  when  the  equilibrium  pressure  p  increases, 
the  same  sensitivity  cannot  be  achieved  thoughout  the  radius  range  of  interest 
with  the  same  diffuser  geometry  which  can  be  analyzed  with  the  methods 
contained  herein. 

Saturation  and  Throat  Section  Volume 

As  pointed  out  briefly  in  the  introduction,  the  problem  of  saturation 
poses  a  significant  limiation  to  the  operation  of  C.S.M. 's.  It  also 
introduces  a  systematic  error,  since  the  unsteady  pressure  perturbations  do  to 
the  presence  of  a  growing  or  collapsing  bubble  in  the  cavitation  region  of  the 
C.S.M.  can  prevent  the  normal  growth  of  other  neighbouring  nuclei. 
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To  control  such  an  error  the  probability  of  finding  more  than  one 
unstable  nucleus  in  the  low  pressure  volume,  of  the  C.S.M.  must  be 

minimized.  When  the  nuclei  concentration  is  uniform,  the  occurrence  of  m 
unstable  nuclei  in  V^p  is  expected  to  follow  the  Poisson  distribution: 


Pr(m) 


m! 


(3) 


where  *  n(RQ)  V^p  is  the  average  number  of  unstable  nuclei  in  V^p.  There¬ 
fore  the  cumulative  probability  of  observing  two  or  more  unstable  nuclei  in 

(bubble  interference  probability,  see  Fig.  6)  is: 

-N  . 

P  -  1  -  Pr (0)  -  Pr(l)  *  -  (1  +  Nab)  e  (4) 

It  is  now  possible  to  express  the  low  pressure  volume  of  the  C.S.M.  as  a  func¬ 
tion  of  the  equilibrium  bubble  radius  and  of  the  bubble  interference  probabil¬ 
ity.  The  results  of  Fig.  7  show  that,  in  the  radius  range  10  to  100  pm,  if 
the  error  is  to  be  reduced  to  an  acceptable  level,  the  volume  of  the  cavita- 
tion  region  of  the  C.S.M.  should  not  exceed  about  50  mm  . 


Bubble  Dynamic  Response  and  Throat  Section  Geometry 

The  problem  of  the  response  of  an  isolated  bubble  in  the  throat  of  the 
C.S.M.  is  now  addressed.  The  solution  of  such  a  problem  in  terms  of  the 
detection  time,  T^^#  md  detection  length,  L(jet  (respectively  defined  as  the 
time  and  the  length  necessary  for  the  bubble  to  grow  to  a  detectable  radius 

s*z®»  ^det^'  b®  as®d  to  estimate  the  minimum  length  of  the  low  pressure 

volume,  thus  completing  the  geometrical  definition  of  this  component. 
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A  number  of  assumptions  are  made  in  order  to  construct  a  simplified  set 
of  equations  which  nevertheless  models  the  interactions  between  the  bubble  and 
the  liquid.  First  the  flow  is  supposed  to  be  one~dimensional ,  ideal, 
incompressible  and  the  relative  motion  of  the  bubble  with  respect  to  the 
liquid  is  neglected.  The  bubble,  moving  along  the  centerline  of  the  duct  with 
instantaneous  position  x^(t),  remains  spherical  and  its  radius  R(t)  is 
determined  by  the  Rayleigh-Plesset  equation  [7]: 


P„  " 


Pb  =  P 


3 

2 


(5) 


where  S  is  the  surface  tension  and  diffusive  and  thermal  effects  are  neglected 
(py  and  Pg  are  constant).  Here  p^  is  the  external  pressure  driving  the  bubble 
volume  changes  and  it  is  assumed  to  be  about  equal  to  the  pressure  of  the 
liquid  in  the  proximity  of  the  bubble.  Then,  from  the  continuity  and  unsteady 
Bernoulli's  equations  for  the  flow  upstream  (subscript  u)  and  downstream 
(subscript  d)  the  bubble: 


Qu  "  %  +  2  '«  R2  S  "  0 


Qd  -  Qb  -  2  *  K2  jf  -  0 


u  dt 


2  \2 


dQj 


2  J2 


x  -51  +  flLLf*  +  i  .0 

d  dt  p  2  2.2 

*b  Ad 


where  _  A(x^),  pQ  and  pex  are  constant  and: 


(6) 

(7) 

(8) 

(9) 
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dx 
A(x ) 


dx 
A  (x ) 


(10) 


I 

a 


x6 


*b 


Equations  (6)  through  (9)  are  then  transformed  by  introducing  perturba¬ 
tion  quantities  with  respect  to  the  steady  state  solution  and  linearized  for 
small  changes  in  the  volume  fluxes.  Finally,  if  I  ~  /  A^>  = 

ut  t  /  At  <  <  Iu  +  1^  ,  which  is  essentially  equivalent  to  consider  the  dynam¬ 
ics  of  a  bubble  entering  a  semi-infinite  pipe  of  constant  cross-sectional  area 

A^,  the  perturbation  and  Rayleigh-Plesset  equations  can  be  reduced  to  the 
following  2nd  order  differential  equation  for  the  bubble  radius  R(t): 


2  n  p  ut 


2  t  R 


2  dh 

dt2 


+  4  t  R 


(ft) 


*ft 


+  p  R 


dt2 


i-  (ft) 


(11) 


M  s 

Pv  -  Pg  \Tj  +  2  R  +  pt 


where  pt  and  u^  are  the  steady  state  pressure  and  velocity  at  the  throat.  The 
detection  time,  T^,  and  detection  length,  computed  ,o-  a  detection 

radius  R£jet  -  150  pm.  and  a  throat  diameter  Dt  ■  1  mm  by  numerically  integrat¬ 
ing  the  above  equation,  are  shown  in  Fig.  8  and  Fig.  9.  The  maximum  in  the 
results  separates  two  regimes:  in  the  lower  radius  range  the  inertial  effects 
of  the  displaced  liquid  dominate,  while  in  the  upper  one  the  effect  of  bubble 
gas  content  prevails,  therefore  accelerating  the  bubble  dynamic  response. 

Note  that,  in  the  radius  range  of  interest,  the  detection  length  is  almost 
independent  of  the  equilibrium  bubble  pressure  and  that  the  length  of  the 
throat  section  should  not  be  shorter  than  5  to  10  mm.  Also  note  that  the 
solution  of  the  above  equation  depends  on  A^f  the  throat  cross-sectional  area. 
Therefore,  if  the  low  pressure  volume  is  initially  fixed  in  order  to  properly 
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limit  the  bubble  interaction  probability,  some  iteration  on  the  throat  diame¬ 
ter  may  be  required  to  completely  define  the  geometry  of  the  C.S.M.  low  pres¬ 
sure  volume. 

Cavitation  Rate  and  Unstable  Nuclei  Concentration  Measurement 

Once  the  throat  cross-sectional  area  has  been  determined,  the  volume 
flux:  q  =  u  A.  and  the  mean  cavitation  event  rate  =  Q  n(R  )  can  be 

t  t  CO 

computed.  Results  shown  in  Fig.  10  and  Fig.  11  for  a  throat  diameter  Dt  =  1  mm 
indicate  that,  in  the  radius  range  from  10  to  100  pm,  as  many  as  a  few  hun¬ 
dreds  of  cavitation  events  per  second  are  to  be  expected. 

For  a  given  pressure  at  the  throat,  the  corresponding  concentration  of 
unstable  nuclei  can  be  ertimated  from  the  count  of  the  cavitation  bubbles 
observed  during  the  sampling  time.  In  steady  state  conditions,  if  the  liquid 
quality  of  the  sample  is  constant,  the  occurrence  of  m  cavitation  events  in 
the  time  interval  (sampling  time)  is  expected  to  follow  the  Poisson 
distribution  with  parameter  ?  =  \)c  Ts,  equal  to  the  average  number  of  events 
in  the  same  time.  Since  can  be  conveniently  measured,  the  problem  of 
extimating  \)^  now  reduces  to  the  estimation  of  £. 

It  is  known  from  the  theory  of  probability  distributions,  [11],  that  the 
best  estimator  of  the  Poisson  parameter  is  the  observed  frequnecy  m,  whose 
variance  is:  V[m]  *  £  *nd  (1  —  u)  confidence  limits  are  given  by: 

h  ‘  2  Xl/2l2‘)  1  h  '2ia/2<2mt2)  (12> 

where  Xp(q)  is  the  p-percentile  of  the  X2  distribution  with  q  degrees  of  free¬ 
dom  (see  Fig.  12).  Thus,  one  way  to  assess  the  relative  indetermination  in  the 
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estimate  of  £  consists  in  normalizing  the  confidence  interval  with  respect  to 

the  corresponding  observed  frequency  (see  Fig.  13,  where  a  =  .20).  Another 

less  rigorous  method  uses  the  estimated  root  mean  square  (r.m.s.)  error  of  the 

observed  frequency  \|  m  ~  \  |  V[m]  instead  of  the  confidence  interval.  The 

results  of  this  second  case  are  presented  in  Fig.  14  and  comparison  with  Fig.. 

13  shows  that  the  two  methods  are  essentially  equivalent  for  large  observed 

frequencies  (m  >  30) ,  apart  from  a  scaling  factor.  Indeed  the  agreement  could 

be  made  closer  by  simply  adjusting  the  confidence  level  1  -  a.  Because  of  its 

simplicity,  the  second  method  is  used  here.  Then,  if  e  is  the  estimated  r.m.s. 

relative  error  in  the  measurement  of  \)  (or  any  other  linearly  related  quan- 

c 

2 

tity,  like  n(RQ)  and  N(RQ)),  the  necessary  sampling  time:  T#  =  1  /  \)c  6 

liquid  volume:  V  =  Q  T  are  plotted  in  Fig.  15,  16  and  17.  From  them  it 
s  s 

appears  that  the  liquid  sample  volume  and  sampling  time  required  for  one  sin¬ 
gle  measurement  of  unstable  nuclei  concentration  of  acceptable  precision  are 
at  most  of  the  order  of  10  dm  and  10  sec  respectively.  Of  course,  the 
overall  sampling  time  and  volume  of  the  liquid  necessary  for  the  evaluation  of 
the  nuclei  number  distribution  will  depend  on  the  total  choice  of  number  of 
points  where  the  measurement  of  the  unstable  nuclei  concentration  is  carried 
out . 

Viscous  Effects  and  Separation 

It  was  previously  mentioned  that  one  of  the  most  difficult  problems  of 
C.S.M. 's  consists  in  measuring  the  throat  pressure  with  the  required  preci¬ 
sion.  Since  direct  measurements  are  virtually  impossible,  indirect  methods 
must  be  used.  Most  naturally  the  throat  pressure  may  be  deduced  from  velocity 
measurements  by  using  Bernoulli's  equation.  An  unavoidable  problem  associated 
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with  this  technique  is  that  the  throat  pressure,  being  inherently  small 
compared  to  the  total  and  kinetic  pressures,  is  expressed  as  the  difference  of 
two  almost  equal  quantities.  Thus,  relative  errors  in  the  evaluation  of  such 
quantities  lead  to  a  larger  relative  error  for  the  throat  pressure. 

In  existing  C.S.M.  's,  [1],  [2],  [3],  [4],  the  velocity  measurement  was 
carried  out  in  the  inlet  section  of  the  duct,  upstream  of  the  contraction. 

This  method  appears  to  be  rather  unsatisfactory  for  various  reasons.  First, 
any  intrinsic  error  in  the  measurement  of  the  upstream  velocity  or  due  to 
inaccurate  estimate  of  boundary  layer  effects  is  amplified  proportionally  to 
the  contraction  ratio  when  the  throat  velocity  is  derived  from  continuity 
arguments.  This  limitation  is  particularly  severe  since  in  C.S.M. 's  the 
contraction  ratio  should  be  as  large  as  possible  in  order  not  to  appreciably 
affect  the  condition  of  the  sampled  liquid  during  its  flow  in  the  inlet  sec¬ 
tion.  Morover,  it  is  generally  necessary  to  make  the  inlet  section  longer, 
which  also  is  likely  to  affect  the  conditions  of  the  sampled  liquid,  espe¬ 
cially  when  a  second  venturi  of  smaller  contraction  ratio  is  used  to  monitor 
the  velocity. 

These  difficulties  are  overcome  if  the  velocity  measurement  is  carried 
out  at  the  throat.  In  consideration  of  the  various  requirements  and  limita¬ 
tions  involved,  the  use  of  a  laser  Doppler  velocimeter  (L.D.V. )  presently 
seems  to  be  the  best  solution.  But,  for  the  pressure  to  be  deduced  from 
velocity  measurements,  the  presence  of  a  potential  core  at  the  C.S.M.  throat 
must  be  demonstrated.  Viscous  effects  have  therefore  been  estimated  by  match¬ 
ing  an  axisymmetric  boundary  layer,  [10],  with  the  quasi  one-dimensional  ideal 
flow  solution  based  on  the  duct  geometry,  [8],  [9],  The  duct  radius  B(x)  has 
been  assumed  to  be  continuous  with  continuous  first  and  second  derivatives,  as 
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required  by  fluid  mechanical  considerations  to  avoid  separation.  In  the 
contraction  and  diffuser  inlet  sections  B(x)  was  supposed  to  have  sinusoidal 
second  derivative.  The  inlet  and  throat  sections  were  taken  to  be  cylindrical 
and  the  diffuser  exit  section  to  be  conical,  of  given  semi-aperture  angle. 

The  length  of  each  section  has  also  been  specified  in  terms  of  the  inlet  duct 
radius. 


The  results  of  such  boundary  layer  calculations  obtained  for  a  prototype 
venturi  made  of  acrylic  plastic  with  Reynolds  number  based  on  the  inlet  radius 
and  average  velocity  Re^  *  700  (the  lowest  expected  value  during  operation) 
are  summarized  in  Fig.  18  through  22.  The  geometrical  characteristics  of  the 
above  prototype  venturi  are:  L  =  11.4  nun,  Lc  =  10.4  mm,  Lt  -  7.7  mm,  Ldi  =  3 
=  6.6  mm,  0d  =  .80  deg,  Dt  *  1  mm,  Cf  =  100,  DCI  /  Dt  =  1.23  and  the 
duct  shape  profile  is  shown  in  Fig.  18,  where  the  scales  on  the  two  axes  are 
different.  The  boundary  layer  displacement  and  momentum  thicknesses  plotted 
in  Fig.  19  (normalized  with  respect  to  the  inlet  radius)  and  in  Fig.  20 
(normalized  with  respect  to  the  local  duct  radius)  clearly  demonstrate  the 
presence  of  a  potential  core  in  the  throat  section  of  the  C.S.M.  In  Fig.  21 
the  pressure  coefficients  at  the  centerline,  C  and  at  the  boundary  layer 
edge,  Cped>  are  shown  as  functions  of  the  axial  coordinate  of  the  duct.  Here 
the  pressure  is  normalized  with  respect  to  the  inlet  pressure  and  velocity. 
Note  the  very  steep  pressure  drop  in  the  contraction  and  the  milder  recovery 
in  the  diffuser. 


Two  kinds  of  problems  may  be  expected  in  the  venturi  as  a  consequence  of 

•  • 

viscous  effects.  If  the  contraction  is  too  steep,  Tay lor-Go rtler  instability 
may  develop,  and,  if  too  large  negative  pressure  gradients  are  present,  flow 
separation  may  occur  at  the  contraction  inlet,  at  the  throat  inlet  or  in  the 


323 


diffuser.  Preliminary  tests  on  the  prototype  venturi  of  Fig.  18  showed  no 
evidence  of  any  of  the  above  mentioned  phenomena,  with  the  exception  of 
intermittent  rotating  separation  in  the  diffuser,  a  problem  reported  by 
Oldenziel  as  well  [1].  In  view  of  this,  the  boundary  layer  correlation  param- 
eter  X  =  0  U'/  \)  (where  0  is  the  momentum  thickness,  0'  is  the  streamwise 
derivative  of  the  boundary  layer  edge  velocity  and  is  the  kinematic  viscos¬ 
ity)  and  the  Stratford  laminar  boundary  layer  separation  criterion,  S  [12] 

c 

have  been  plotted  in  Fig.  22.  Separation  is  predicted  when  X  =  -.09  or 
Sc  -  0.  The  two  methods  agree  well  with  each  other  and  also  with  the  observed 
location  of  separation  in  the  prototype  venturi,  thus  indicating  that  less 
steep  diffusers  should  be  used. 

SUMMARY  AND  CONCLUSIONS 

The  operation  of  C.S.M. 's  has  been  analyzed  in  order  to  investigate  the 
possibility  of  extending  their  use  to  determine  the  concentration  of  unstable 
cavitation  nuclei  in  a  liquid  as  a  function  of  the  applied  tension,  from  which 
the  nuclei  number  distribution  can  be  deduced. 

By  means  of  various  assumptions  and  approximations,  a  simplified  model 
has  been  derived  which  adequately  represents  the  operation  of  C. S.M. 's  in 
terms  of  a  limited  number  of  relevant  parameters.  The  investigation  of  their 
mutual  relations  led  to  a  functional  description  of  C.S.M.  operation  which 
enables  one  to  address  the  problem  of  their  design  in  a  systematic  and  organic 
way. 


The  above  model  has  been  applied  to  simulate  C.S.M.  operation  for  the 
case  of  a  typical  nuclei  number  distribution  in  water.  The  results  can  be 
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summarized  as  follows: 


a)  the  pressure  in  the  C.S.M.  cavitation  region  is  an  extremely  sensitive 
parameter  and  cannot  be  measured  directly  with  the  necessary  accuracy; 

b)  the  indirect  measurement  of  the  pressure  in  the  cavitation  region  is 
best  accomplished  by  the  measurement  of  the  potential  core  velocity  at 
the  C.S.M.  throat  by  means  of  a  L.D.V.; 

c)  the  existence  of  a  potential  core  at  the  throat  has  been  demonstrated 
by  estimating  boundary  layer  effects; 

d)  for  the  specific  nuclei  number  distribution  considered,  the  errors  and 
limitations  due  to  saturation  can  be  reduced  by  appropriately  decreasing 
the  volume  of  the  cavitaticn  region; 

e)  the  time  response  of  the  cavitating  nuclei  imposes  a  lower  bound  for 
the  length  of  the  C.S.M.  throat  section; 

f)  the  flow  can  sustain  without  separation  or  fluid  mechanical 
instabilitites  the  very  large  and  steep  contraction  suggested  by  design 
considerations,  but  is  extremely  sensitive  to  separation  in  the  diffuser; 

g)  in  the  case  of  the  nuclei  number  distribution  used  as  an  example,  it 
appears  that  C.S.M. 's  can  be  designed  which  are  capable  of  detecting  gas 
bubbles  having  an  equilibrium  radius  larger  than  about  10  |un  within  the 
errors  typical  of  other  alternative  techniques  currently  available; 

h)  due  to  the  considerable  complexity  of  the  phenomena  involved,  the 
interpretation  of  the  data  obtained  from  C.S.M. 's  requires  a  careful 
analysis. 


325 


ACKNOWLEDGEMENT 


The  authors  would  like  to  thank  Carol  Kirsch  for  her  help  in  organizing 
this  manuscript.  This  work  was  supported  by  Naval  Sea  Systems  Command  General 
Hydromechanics  Research  Program  Administered  by  the  David  Taylor  Naval  Ship 
Research  and  Development  Center  under  Contract  No.  N00014-75-C-0378. 


326 


REFERENCES 


1.  Oldenziel,  D.  M. ,  "A  new  instrument  in  cavitation  research:  the  cavita¬ 
tion  susceptibility  meter",  ASME  Journal  of  Fluids  Engineering,  104,  June 
1982,  pp.  136-142. 

2.  Oldenziel,  D.  M. ,  Jansen,  R.  H.  J. ,  Keller,  A.  P.  ,  Lecoffre,  Y. ,  and  van 
Renesse,  R.  L. ,  "Comparison  of  instruments  for  detection  of  particles  and 
bubbles  in  water  during  cavitation  studies",  Proceedings  of  Symposium  on 
Operating  Froblems  of  Pump  Stations  and  Power  Plants,  IAHR,  .Amsterdam, 
September  1982. 

3.  Oldenziel,  D.  M. ,  "Utility  of  available  instruments  during  cavitation 
tests".  Proceedings  of  Symposium  on  Operating  Problems  of  Pump  Stations 
and  Power  Plants,  IAHR,  Amsterdam,  September  1982. 

4.  Godefroy,  H.  W.  H.  E. ,  Jansen,  R.  H.  J.  ,  Keller,  A.  P.  ,  and  van  Renesse, 
R.  L. ,  "Comparison  of  measuring  and  control  methods  of  the  water  quality 
with  respect  to  cavitation  behaviour".  Delft  Hydraulics  Laboratory  Publi¬ 
cation.  January  1981. 

5.  Knapp,  R.  T.  ,  Daly,  J.  W.  and  Hammit,  F.  G. ,  "Cavitation",  McGraw  Hill, 
1970. 

6.  van  Wijngaarden,  L.  ,  1972,  "One-dimensional  flow  of  liquids  containing 
small  gas  bubbles",  Ann.  Rev.  Fluid  Mech. ,  Vol .  4,  pp. 3 69-3 96. 

7.  Plessot,  M.S,  and  Prosperetti,  A.,  "Bubble  dynamics  and  cavitation",  Ann. 
Rev.  Fluid  Mech.,  1977,  Vol.  9,  145-85. 


327 


8.  Chmielewski,  G.  E.  ,  "Boundary  layer  consideration  in  the  design  of 
aerodynamic  contractions",  J.  of  Aircraft,  Vol .  11,  No.  8,  August  197  4, 
pp.  43  5-438. 

9.  Mikhail,  M.  N. ,  "Optimum  design  of  wind  tunnel  contractions",  AIAA  Jour¬ 
nal,  Vol.  17,  No.  5.  May  197  8,  pp.  471-477. 

10.  White,  F.  M.  ,  "Viscous  fluid  flow".  Me  Graw  Hill,  197  4. 

11.  Brownlee,  K.  A.,  "Statistical  theory  and  methodology  in  science  and 
engineering",  John  Wiley  4  Sons,  1960. 

12.  Stratford,  B.  S. ,  "Flow  in  laminar  boundary  layer  near  separation". 
Aeronautical  Research  Council,  RAM  3002,  November  1954. 

13.  Le  Goff,  J.  P.  and  Lecoffre,  Y. ,  "Nuclei  and  cavitation",  14th. 

Symposium  on  Naval  Hydrodynamics,  National  Academy  Press,  1983  ,  pp.  215- 
242. 


328 


c n 

3 

< 

X 

X 

UJ 


x  3 
LU  < 

D  5 
£  2 
U.  O 


b 


329 


igure  1.  Axi symmetric  duct  geometry  and  nomenclature 
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Figure  4.  Throat  velocity,  u^  [OT],  v/ a  babble  equilibrium  radius,  R  [RO] 
at  various  equilibrium  pressures:  p  =1  bar  (lover  curve),  p  =10  bar 
(intermediate  curve)  and  pQ  =  30  bar  (upper  curve)  for  S  =  .0^3  N/m  (surfac 
tension),  p  =  1919  Pa  (vapor  pressure)  and  p  =  1000  kg/m^  (liquid  density) 


1E  +  00 


333 


1E  +  00 


co 


33^ 


Figure  6.  Bubble  interference  probab  iity,  0  [BETA],  i.e  probability  of  find¬ 
ing  more  than  one  bubble  in  the  cavitation  region  of  the  duct,  v/s  the  averagi 
number  of  unstable  bubbles  in  the  same  region,  N  .  [NUB],  Such  probabilities 
are  calculated  assuming  that  the  occurrence  of  cavitation  nuclei  in  a  given 
volume  is  represented  by  a  Poisson  process. 
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Figure  10.  Volume  flow  rate.  Q,  v/ s  equilibrium  bubble  radius,  R  [R0]  ,  at 
various  equilibrium  pressures:  p  =1  bar  (lower  curve),  p  =  10  §ar 
(intermediate  curve)  and  pQ  =  30°bar  (upper  curve)  for  S  =°.073  N/m  (surface 
tension),  p  =  1919  Pa  (vapor  pressure),  p  =  1000  kg/m  (liquid  denisty)  and 
D  =  1  mm  (throat  diameter). 
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Figure  12.  Poisson's  distribution  parameter  estimator,  5  [CSIEST],  and 
confidence  limits,  [CSI1]  and  L  [CSI2] ,  v/ s  sample  statistical  frequence 
m  [M] .  Here  the  confidence  level  Is:  1  -  a  =  .8. 
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Figure  13.  Poisson's  distribution  parameter  relative  confidence  interval 
($2  ~  5j)/b  [RCI],  v/ s  sample  atatistical  frequency,  m  [M] .  Here  the 
confidence  level  is:  1  —  a  =  .8  and  the  confidence  interval  is  normalized 
w.r.t.  the  sample  frequency. 
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Figure  14.  Poisson’s  distribution  parameter  estimated  relative  mean  quadratic 
error,  s  [ERMQE] ,  v/s  sample  statistical  frequency  m  [ M ]  .  Here  the  best  esti¬ 
mator  of  the  Poisson's  distribution  parameter  (the  sample  frequency)  i«  used 
to  both  evaluate  the  r.m.s.  error  and  to  normalize  it. 


343 


Figure  15.  Sanpling  time,  [TSAM] ,  v/ s  average  cavitation  event  rate, 
[NUC] ,  for  various  values  of  the  relative  estiaated  r.m.s.  error:  e  =  0.1 
(upper  curve),  .05  (intermediate  curve)  and  .20  (lower  curve). 
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LIQUID  SAMPLE  VOLUME,  VSAM,  V/S  EQUILIBRIUM  BUBBLE  RADIUS.  R0 
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calculations  are  based  on  a  quasi  one-dinensional  node!  of  the  flow  in  the 
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